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This paper deals with the study of the accretion of a generalized Chaplygin gas with equation 
of state p = —AIp“‘ onto wormholes. We have obtained that when dominant energy condition is 
violated the size of wormhole increases with the scale factor up to a given plateau. On the regime 
where the dominant energy condition is satisfied our model predicts a steady decreasing of the 
wormhole size as generalized Chaplygin gas is accreted. Our main conclusion is that the big trip 
mechanism is prevented in a large region of the physical parameters of the used model. 
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Several astronomical and cosmological observations, ranging from distant supernovae la Q to the cosmic microwave 
background anisotropy Q , indicate that the universe is currently undergoing an accelerating stage. It is assumed that 
this acceleration is due to some unknown stuff usually dubbed dark energy, with a positive energy density p > 0 and 
with negative pressure p < —{l/3)p. 

There are several candidate models for describing dark energy, being the cosmological constant. A, by far, the 
simplest and most popular candidate Other interesting models are based on considering a perfect fluid with given 
equation of state like in quintessence^^ , K-essence 0 or generalized Chaplygin gas Q, models. Note 

that there are also other candidates for dark energy based on brane-world scenarios [I^l and modified 4-dimensional 
Einstein-Hilbert actions M , where a late time acceleration of the universe may be achieved, too. 

One of the peculiar properties of the resulting cosmological models is the possibility of occurrence of a cosmic 
doomsday, also dubbed big rip H. The big rip appears in models where dark energy particularizes as the so-called 
phantom energy for which the dominant energy condition is violated, so that p+p < 0. In these models the scale factor 
blows up in a finite time because its cosmic acceleration is even larger than that induced by a positive cosmological 
constant. In these models every component of the universe goes beyond the horizon of all other universe components 
in finite cosmic time. It should be noted that the condition p -|- p < 0 is not enough for the occurrence of a big rip, 
i.e., if one consider an universe filled with phantom generalized Chaplygin gas, one can avoid the big rip M (see also 
hd). Other peculiar properties of phantom energy are that it can make the exotic substance that fuel wormholes 
|17|. triggering the possibility of occurrence of a big trip , i.e. if there is a wormhole in an universe filled 

with phantom energy, due to processes of phantom energy accretion onto the wormhole, the size of this wormhole 
increases in such way that the wormhole can engulf the universe itself before it reaches the big rip singularity, at 
least relative to an asymptotic observer. Then, the following question arises: Does an universe filled with generalized 
phantom Chaplygin gas which avoids the big rip singularity also escape from the big trip? We have found in the 
present paper that for a large range of the Chaplygin parameters no big trip is predicted by our model, though there 
still exists sufficient room in the physically allowed parameter space to not exclude the possibility for the occurrence of 
such a strange causality disruption. Thus, generalized Chaplygin gas can be seen to have several interesting potential 
properties, as it avoids big rip it may be used as the stuff to construct wormholes |2l|. it prevents the universe 
to be engulfed by a black hole |^ . and, finally it may also circumvent the big trip problem. 

We starts by first reviewing the accretion formalism first considered by Babichev, Dokuchaev and Eroshenko 
(see also H^), generalizing it to the case of wormholes. Throughout this paper we shall use natural units so that 
G = c = 1. The Morris-Thorme static space-time metric of one wormhole is given by 

ds^ = —-I- + r^{d9^ -|- sin^ Odcf^), (1) 
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where $(r) is the shift function and K{r) is the shape function. We model the dark energy in the wormhole by the 
test perfect fluid with a negative pressure and an arbitrary equation of state p{p), with the energy-momentum tensor 


Tfii, = {p + p) u^Uv + pgfj,„, (2) 

where p is the pressure, p is the energy density, and = dx^/ds is the 4-velocity with = —1. The zeroth (time) 

component of the energy-momentum conservation law = 0 can then generally be written as 





(P + P) 


dt dr 
ds ds 


+ (p + p) 


$'(r) 


K'{r)r — K{r) 

2r2 (l - j 


2 dt dr 
r ds ds 


( 3 ) 


This expression should now be integrated. The integration of Eq. m gives then. 


uM (p _|_ p) ^y2 I _ = c, (4) 

where u = dr/ds, and M is the exotic mass of the wormhole which, following the procedure of Ref. |^, has been 
introduced to render the integration constant C to have the dimensions of an energy density (note that we are using 
natural units), and, without any loss of generality for our present purposes, we have adhered to the case where = 0. 

Another integral of motion can be derived by using the projection of the conservation law for energy-momentum 
tensor along the four-velocity, i.e. the flow equation = 0. For a perfect fluid, this equation reduces to 


u^P,ri + {P + P)u‘^-,r- = 0 - 

The integration of Eq. gives the second integral of motion that we shall use in what follows 
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P+P = —A, 


( 6 ) 


where rt < 0 in the case of a fluid flow directed toward the wormhole, and A is a positive dimensionless constant. 
Eq. (O gives us the energy flow induced in the accretion process. From Eqs. 0) and 10 one can easily get 


(p + p) 



K(r) _ fp 
- e p°° 


dp 

P+P 


= C2 


( 7 ) 


where C 2 = —CfA = A {p{poo) + Poo), with A a positive constant. 

The rate of change of the exotic mass of wormhole due to accretion of dark energy can be derived by integrating 
over the surface area the density of momentum Tq’’, that is 


M = - 


/ 


To^dS, 


with dS = r^ sin 9d9d(j),. Using Eqs. 0, ® and 0 this can be rewritten as |2fl| 


( 8 ) 


M = -AttDM'^ 



K{r) 


\p {Poo) 4” Poo] , 


(9) 


with the constant D = A A > 0. For the relevant asymptotic regime r ^ 00 where the big trip occurs, the rate M 
reduces to 


M = {p + p) . 


( 10 ) 


We see then that the rate for the wormhole exotic mass due to accretion of dark energy becomes exactly the negative 
to the similar rate in the case of a Schwarzschild black hole, asymptotically. For current quintessence models, the use 
of a scale factor i? = i?o (l + V67J'Po(l + w){t — for li; < —1 (which corresponds to the solution of the 
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general equation —3H{1 + w) = 2H/H^ with ' = d/dt and H = R/R for r(;=constant ^3), gives arise to a singular 
behaviour of the size of the wormhole before reaching the big rip, that is the big trip. 

We shall derive now the expression for the rate M in the case of a generalized Chaplygin gas. This can be described 
as a perfect fluid with the equation of state 


p = -Ach/p°', ( 11 ) 

where Ach is a positive constant and a > —1 is a parameter. In the particular case a = 1, the equation of state uu 
corresponds to a Chaplygin gas. The conservation of the energy-momentum tensor implies 


P = 







( 12 ) 


with B = — Ach) R = R{t) is the scale factor and the subscript “0” means initial value. Now, from 

the Friedmann equation we can get 


R = 




i?3(l+“) ) 


l/[ 2 (l-ra)] 


Hence, from Eqs. (tTUli-itnii we obtain 


(13) 


M = 


Mo 


1-DMo 




(14) 


For the case where the dominant energy condition is preserved, i.e. i? > 0, we obtain that M decreases with time and 
tends to a constant value. On the other hand, M is seen to decrease more rapidly as parameter a is made smaller. If 
the dominant energy condition is assumed to be violated, i.e. H < 0, as phantom energy is assumed to require m, 
then M increases with time, with M tending to maximum, nonzero constant values. Making \B\ or a smaller, makes 
the evolution quicker. 

When time goes to infinity, then the exotic mass of wormhole approaches to 


M = 


Mo 


l-DMo 





(15) 


that is a generally finite value both for B > 0 and B < 0. Thus, at first sight it could be thought that, unlike 
what happens in phantom quintessence models, the presence of a generalized Chaplygin gas precludes the eventual 
occurrence of the big trip phenomenon. However, such a conclusion cannot be guaranteed as the size of the wormhole 
throat could still exceed the size of the universe during its previous evolution. Note that for e.g. wormholes with 
zero tidal force, one can consider that the exotic matter is confined into an arbitrarily small region around the 
wormhole throat, then the radius of wormhole throat becomes roughly proportional to its mass. That is to say, the 
question remains whether the wormhole would grow eventually rapidly enough or not to engulf the universe during the 
evolution to its final classically stationary state. That question should be settled down before reaching a conclusion on 
the possibility of the big trip in these models. Actually, in order for avoiding a big trip, the following two conditions 
are also required: (i) i? M along the entire evolution, and (ii) that N{R) = R/M be always an increasing function 
along that evolution. The first of these conditions implies that the function 


f{M) = M - MMoD\l — 
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be nonvanishing everywhere. We analyze this question by taking the zeros of the second derivative 


/"(M) = 


d^f{M) 

dM^ 


= V^DB 


Mo 


M3(l+a) + l 


A. 


B 


ch 


M3(1+“) 


-2q!-i r 
2(l + «) 


1 — 3 (1 -l- Of) -l- 


W{2a+1) 

2M3(i+“) 


B 


^ch 


M3(1+“) 


-1 


(17) 





















Now, by taking into account that physically, 
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^ \ 3(l + c) 

M > Mo > ( 


(18) 


The second inequality meaning that, during its evolution, the radio of the universe should be larger than its initial 
size. Eq. 03 can be reduced to imply for the zeros 


1 — 3 (1 + a) + 


3B(2a + l) 
2M3(i+“) 


A. 
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ch 




-1 


= 0 , 


(19) 


whose solution would read = —i?/[2(2 + 3a)Ac.h]- Nevertheless, from condition m we have that 0 < 

2(2 + 3a) < 1 which, in turns, implies that there could be at most three crossing points, as this interval allows for 
values of a within the range of the generalized Chaplygin gas models, i.e. a > — 1. For a-values outside the interval 
0 < 2(2 + 3a) < 1, but still inside the Chaplygin range, one might expect just two points at most. It follows from 
this analysis that a big trip could in principle happens. 

As to condition (ii), we obtain from Eqs. (Cnj and m that 


N{R) 


_ - _i?3(l+a) + l ^ 

BM^Dy/^ 


1 - MoD 




( 20 ) 


From Eq. 1201) it follows that N{Ro) > 1 because the size of the wormhole must be quite smaller than that of the 
universe initially. Thus, an always increasing N{R) should necessarily imply that the scale factor increased more 
rapidly than M did, so preventing any big trip to occur. By differentiating N{R) with respect to R and taking into 
account that i? < 0 for Chaplygin phantom, it can be checked that N{R)/dR > 0 in the general case that a does not 
reach values sufficiently close to —1; that is, it follows from Eq. 03 that inside the interval 
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( 21 ) 


a big trip would still take place. 

On the other hand, the question still remains on what happens with the grown-up wormhole once it has reached its 
maximum, final size. Since the wormhole size tends to become constant at the final stages of its evolution and it is 
rather a macroscopic object, it would be subjected to chronology protection 123 ■ In fact, one expected that vacuum 
polarization created particles which catastrophically accumulated on the chronology horizon of the wormhole making 
the corresponding renormalized stress-energy tensor to diverge and hence the wormhole would disappear. 

It would be quite interesting to probe the region of parameter space (a, Ach, Hq, fix, fl^) allowed by current obser¬ 
vations in order to determine whether there exist any allowed sections leading to a big trip. However, all available 
analysis |2^ |2^ 123,0 are restricted to the physical region where no dominant energy condition is violated. There¬ 
fore, the section described by the interval implied by Eq. necessarily is outside the analyzed regions. One had 
to extend the investigated domains to include values of parameter Ach > 1 to probe the parameters space relative to 
the big trip. In any events the range of a values compatible with a big trip has been seen to be extremely narrow 
and hence the occurrence of the big trip phenomenon in the generalized Chaplygin gas model appears to be highly 
unlikely. 

In this paper we have studied the accretion of a generalized Chaplygin gas onto a wormhole. First, we have reviewed 
the accretion formalism originally considered by Babichev, Dokuchaev and Eroshenko 123 for the case of a wormhole 
1^. We have then applied such a formalism to the generalized Chaplygin gas model. The evolution of exotic mass 
with the accretion of Chaplygin dark energy has been first considered for the case that the dominant energy condition 
is satisfied. It has been seen that in that case the mass decreases with cosmic time. 

If accretion involves Chaplygin phantom energy, then M increases from its initial value, tending to reach a plateau as 
cosmic time goes to infinity. It is obtained that for a wide region of the Chaplygin parameters no big trip is predicted, 
contrary to what happens in quintessence and K-essence dark energy models. However, as far as the Chaplygin regime 
tends to match the quintessence regime, but presumably still within the Chaplygin region, the possibility for a big 
trip at a finite time in the future is not excluded. Finally, we also argued that the fate of the final wormhole is to be 
destabilized by quantum vacuum processes. To conclude, the generalized Chaplygin gas has several very interesting 
features and may circumvent several singularities that appear in the usual quintessence models. Whether or not the 
above features can be taken to imply that Chaplygin gas is a more consistent component than usual quintessential or 
K-essence dark energy component with in < —I is a matter that will depend on both the intrinsic consistency of the 
models and the current observational data and those that can be expected in the future. 
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